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ABSTRACT 

A statistical analysis of texture on the COB frDMR first year sky maps based on the 
genus and spot number is presented. A generalized \ 2 statistic is defined in terms of 
"observable" quantities: the genus and spot density that would be measured by different 
cosmic observers. This strategy together with the use of Monte Carlo simulations of the 
temperature fluctuations, including all the relevant experimental parameters, represent 
the main difference with previous analyses. Based on the genus analysis we find a strong 
anticorrelation between the quadrupole amplitude Qrms-PS and the spectral index n 
of the density fluctuation power spectrum at recombination of the form Qrms-PS = 
22.2 ± 1.7 — (4.7 ± 1.3) x n [iK for fixed n, consistent with previous works. The result 
obtained based on the spot density is consistent with this Qrms-ps(n) relation. In 
addition to the previous results we have determined, using Monte Carlo simulations, 
the minimum uncertainty due to cosmic variance for the determination of the spectral 
index with the genus analysis. This uncertainty is 8n 0.2. 

Key words: Cosmic Microwave Background Radiation - Cosmology 



1 INTRODUCTION 

The study of texture on Cosmic Microwave Background 
(CMB) maps can be used to constrain scenarios of galaxy 
formation as an alternative technique to the temperature 
correlation analysis. Several statistical quantities have been 
proposed, such as number of spots, contour length, genus 
(total curvature of the iso-temperature contours), and to- 
tal area of excursion above a certain level (Sazhin 1985, 
Bond & Efstathiou 1987, Vittorio & Juszkiewicz 1987, Coles 
1988, Martmez-Gonzalez and Sanz 1989, Gott et al. 1990). 
At large angular scales (> 2°) gravitational potential fluc- 
tuations at the cosmic photosphere must leave an imprint 
on the CMB (Sachs & Wolfe 1967), which is manifested as 
temperature fluctuations on the 2D sphere. 

There have already been several statistical analyses of 
the COBE data. However, due to the great impact that 
these data have on cosmology and the difficulties of its in- 
terpretation caused by the experimental complexities and 
the low signal level, it is important to look at these data 
exhaustively. We propose and use a new statistic, based 
on the topological characteristics of sky maps, which is di- 
rectly related with observable quantities. Analysis of the first 
year COBE-DMB. maps based on the genus have recently 
been carried out by Torres (1994a, 1994b) and Smoot et al. 



(1994). From the first paper a value of n = 1.2 ± 0.3 was 
obtained by fitting the coherence angle of the COBE tem- 
perature maps for a fixed quadrupole-normalized amplitude 
of Qrms-ps = 16fiK, assuming a scale-free primordial spec- 
trum of density fluctuations P(k) cx Q 2 rrns _ PS k n . Smoot et 
al. (1994) allow variations of Q and find a relation between 
Q and n based also on genus analysis. In Smoot et al. (1994) 
the covariance matrix of the genus at different temperature 
thresholds and for several smoothing scales is included in 
the x 2 minimization. The main difference of our analysis 
with the previous works lies in the direct comparison of the 
genus and spot density of the COBE data with the values 
obtained in each realization of a model with no additional 
smoothing of the data. 

The previous analyses have been performed by compar- 
ison between the observed genus and the mean values de- 
duced from the cosmological models. Mean values of model 
parameters are not observable quantities and therefore are 
not appropriate to compare with the observed data (we will 
see, however, that in some cases a statistic based in the 
mean value of a quantity can be interpreted as the mean 
value of the statistic based on realizations of that quan- 
tity; see section 3.1 below). In this work we use a different 
statistical analysis of the genus and spot density data that 
takes into account the ensemble of realizations that would 
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be measured by different cosmic observers. This approach 
has been used by Scarameifa and Vittorio (1993) with the 
temperature correiation function as the observabie quantity. 
However, they do not consider aii the characteristics of the 
COBE-DMB. experiment such as gaiactic cut, pixeiization, 
beam smearing, and smoothing, which is essentiai for a pre- 
cise determination of the values of cosmological parameters 
allowed by the data. 

Improvements in instrument sensitivity are very rapidly 
reducing experimental errors, however even in an ideal noise- 
less experiment the measured parametes will have the uncer- 
tainty due to cosmic variance. We have used Monte Carlo 
simulations to study the extent to which cosmic variance 
obscures the information which can be extracted from the 
genus analysis. A similar consideration for the rms temper- 
ature fluctuations on different angular scales results in min- 
imal uncertainties in the determination of n (White, Krauss 
& Silk 1993). 

In section 2 we present how Monte Carlo simulations of 
the COBE-DMB. maps are performed. The characteristics 
of the new statistical method considered in this paper are 
presented in section 3. The minimum range for the spectral 
index n, implied by cosmic variance, is obtained in section 
4. In section 5 we apply the new method to the genus and 
spot density of the COBE-DMB. data and show the results. 
An interpretation of these results in terms of the coherence 
angle is given in section 6. Finally, a summary of the main 
conclusions is given in section 7. 



2 MONTE CARLO SIMULATIONS 

Simulations that take into account instrument noise, sky 
coverage, galactic cut, smearing, pixeiization scheme and 
the DMR beam characteristics were done for a set of cos- 
mological models. Simulated maps of the cosmic signal were 
generated using an expansion on real harmonics for the tem- 
perature (Smoot et al. 1991) for each one of the 6144 DMR 
pixels: 
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where k = \/2 for to ^ and k = 1 for to = 0; P ; m (cos 8) are 
the Associated Legendre polynomials; the bi tTn coefficients 
are real stochastic and Gaussian distributed variables with 
zero mean and model dependent variance (hi m ) given by 
(Abbott & Wise 1984, Bond & Efstathiou 1987) 
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The weights W t for DMR given by Wright et al. (1994) were 
used. For each realization two maps are generated by adding 
to the cosmic signal the noise corresponding to channels A 
and B of COBE 53 GHz. The noise is determined by instru- 
ment sensitivity and the number of observations per pixel. 
A small beam smearing correction is applied in order to take 



into account the motion of the spacecraft during the 1/2 sec- 
ond integration time. The two maps are added to form the 
j(A+B) map, Gaussian smoothed (<r s = 2.9°) and finally 
its genus is calculated. The algorithm used to compute the 
spot number density and genus is described elsewhere (Tor- 
res 1994a; 1994b). 



3 THE NEW STATISTICAL METHOD 

We define the statistic x% associated to the genus and Xn 
associated to the number of spots as follows: 
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G? is the genus for the fc-th realization at threshold level 
i, and 25 equally spaced threshold levels were used (y : 
-3.0, ...,3.0); G? OBE is the genus for the COBE map at 
threshold i. It is important to notice that Gf is the genus 
as measured by a cosmic observer with an instrumental de- 
vice like that aboard COBE and thus the statistic {Xa) k 
so constructed is based on observable quantities. M %] is the 
covariance matrix: 
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The Mij matrix was calculated with the Monte Carlo real- 
izations. A statistic (Xn) related with the spot number is 
defined in a similar manner as {Xa) k by replacing Gf (and 
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for the realization k at threshold i. So each model is now 
represented by a distribution of values {Xa) k an d in order 
to find the model that is closest to the COBE data we have 
to compare distributions. The simplest way to do this is to 
use the mean value of the distribution < Xg > (although 
other choises like the mode can be used they are more sensi- 
tive to the noise due to the limited number of realizations ). 
Notice that considering the mean value, < Xg >> i s equiv- 
alent to a statistic defined by replacing the genus for each 
realization Gf in equation (3) by the expected genus of the 
model, except for the constant value 25. Another possibility 
is to use the Kolmogorov-Smirnov statistic K to compare 
the Xg distributions. That we will do at the end of section 
5 in order to check the results. For the moment we will take 
the simplest alternative < Xg > which also requires less 
CPU time to compute. 



4 IDEAL EXPERIMENT 

In this section we consider an ideal experiment, that includes 
all the COBE-DMB. experiment characteristics but assum- 
ing a noise-less radiometer, to elucidate the main properties 
of the new statistical method presented and to study the 
cosmic variance associated to the cosmological parameters 
as derived from the genus and spot number. The noise of the 
COBE-DMB. radiometers is taken into account when apply- 
ing the method to the real data in next section. 

The no ergodic character of the CMB temperature ran- 
dom field on the Cosmic photosphere implies a limitation in 
the accuracy to determine the statistical properties of the 
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field from a single realization. Even in the ideal case that we 
were able to measure the background temperature fluctua- 
tions in every part of the sphere (no galactic cut) and with 
negligible noise we would only be able to measure the ex- 
pected values of the field within a certain error bar (cosmic 
variance). Uncertainties in the two point correlation function 
and in other higher moments due to cosmic variance have al- 
ready been calculated (Scaramella and Vittorio 1990, Cayon, 
Martinez- Gonzalez and Sanz 1991, White, Krauss and Silk 
1993, Srednicki 1993, Gutierrez de la Cruz et al. 1994). 

We calculate the cosmic variance uncertainty when 
analysing the genus of a hypothetical sky map with the 
characteristics of the COBE-DMB. experiment assuming a 
noise-less radiometer and no galactic contamination. In this 
situation the genus is independent of the amplitude of the 
cosmic signal and the only free parameter for the models is 
n. Thus, 7 Monte Carlo data sets were generated, one for 
each value of n in the range 0.0-2.0. The estimated value of 
n was obtained using the generalized \ 2 method described 
above with each realization from model n = 1 taken as the 
input data (i.e. used as the COBE data in eq. 3). So, from 
the whole n = 1 data set we obtain a distribution of n val- 
ues. The dispersion of this distribution gives the uncertainty 
due to cosmic variance: Sn ~ 0.2. 



5 ANALYSIS OF THE COBB-DMR MAPS 

5.1 COBE-BMR maps 

Only data from the most sensitive radiometers were used 
(i.e. the DMR 53 GHz). Before any analysis was done, the 
released maps were processed as follows: 1. The maps were 
converted to thermodynamic temperature scale (a factor of 
1.0742). 2. A monopole and dipole function, including the 
small quadrupolar component of the Doppler effect, was fit- 
ted to the maps (excluding ±30° from the galactic plane) 
and subtracted. The dipole and quadrupole fit is in agree- 
ment with Smoot et al. (1992), thus providing a check for 
the integrity of the data and the analysis software. 3. Finally 
the sum i(A+B) and difference i(A-B) maps were formed 
and Gaussian smoothed (a s = 2.9°) in order to reduce noise. 

5.2 Analysis 

To find the restrictions on Qrms-ps and n imposed by the 
COBE-DMB. genus and spot density, a grid of Monte Carlo 
data sets were generated with n in the range 0-2 (with a 
step of 0.2) and Qrms-ps between 5 and 33 fiK (with a 
step of 2/tK). For each of the simulated i(A+B) maps , the 
statistics Xg an d Xn were calculated as indicated in section 
3. The number of realizations was set to 400, which proved 
to be a safe number after testing for the convergence of the 
results of the relevant quantities. 

For each model (i.e. pair of values Q rms -ps, n) one can 
build up, using the Monte Carlo data set, the probability 
Pg(xg) °f obtaining a realization with its Xg smaller or 
equal to Xg- Thus, models with a high Pg(xg) f° r small 
Xg lm ply that many observers would measure a texture of 
their maps very similar to that of COBE, and the opposite 
for models with flatter Pg(xg)- The same can be done for 
the probability Pn(x. 2 n) based on the spot number. Two 



Figure 1. Distributions of the x% statistic for two models: n = 
0.4, Qrms — PS = 19 (solid) and n = 0.8, Qrms — PS = 19 (dash). 



examples of the distribution of x% f° r models n = 0.4, 
Q rm s-ps = 19/tK and n = 0.8, Q rms -ps = 19/tK are 
presented in Fig. 1. For both models the histogram has a 
maximum near x% ~ 45 and the two distributions are very 
similar (as we will discuss later, there is not a significant 
statistical difference between them). 

In order to establish a simple criterium to compare 
probabilities, and so to choose the model that best fit the 
data, we have used the mean < Xg > (< X%r >) value of 
the realizations for each model. We then search for the min- 
imum < Xg > (< X% >) m the space of model parameters 
(Qrms-ps, n). The first and main result is the anticorrela- 
tion found between our estimates of n and Q rms -ps which 
can be approximately given as a straight line of the form: 
Q rms _ PS = 22.2 ± 1.7 - (4.7 ± 1.3) x n fiK for fixed n. An 
explanation of this anticorrelation in terms of the coherence 
angle is given in the next section. To obtain the previous 
relation between Q rms -ps and n, we first assign an error 
bar to the value of Q rms -ps with minimum < Xg > f° r 
fixed n by using the generalized % 2 method with each real- 
ization from that Q rms -ps taken as the input data. We then 
fit a straight line to the pairs (Q rms -ps,n) with minimum 
< Xg > considering the corresponding error bars. It is clear 
from this analysis that there is a wide range of values of n 
that fit the data equally well (as can be seen in figure 1). 
However it seems that for the < Xg > (< X% >) statistic 
values of n lower than 1 are favored. We have checked the 
stability of the main result with the number of realizations 
by performing an additional set of 400 simulations for mod- 
els which differ from the minimum < Xg > m 1-5 units (a 
number bigger than the typical error for a sampling of 400 
realizations). It is verified that the (Q rms -ps,n) relation is 
mantained. 

The error bars for the cosmological parameters 
Q rms -PS, n were estimated by using each one of the realiza- 
tions for the model n = 0.4, Q rms -ps = 19/tK (which lies 
in the line of degeneracy) as input data replacing COBE's 
data. The result is a 68% uncertainty of ±0.4 for n and 
for Q r ms-ps , which we consider as typical error bars for n 
and Q r ms-ps- 

Similar results were found for the number of spots N 
but for this quantity the generalized % 2 analysis was less 
sensitive to variations in the model parameters. 
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We now test the robustness of the < % 2 > method to 
discriminate among the various cosmological models by com- 
paring it with a more sophisticated method based on the 
Kolmogorov-Smirnov K statistic. The K statistic is used 
to compare the distribution of values (xg)cobe obtained 
from equation (3) with the distribution (xg)i constructed 
by replacing Gf' OBE in that equation by the genus of the 
realization I of the same model, G\, and calculate the statis- 
tic K l G for the two distributions. Operating in the same way 
with all the realizations we then obtain a distribution of K l G 
values for that model and we use the mean value < Kg > to 
compare different models. With the new < Kg > minima we 
are able to approximately recover the relation (Q rms -ps, n) 
given above. However we do not find the trend found with 
< Xg > favoring n < 1 models. Given the relatively poor 
signal to noise ratio of the DMR first year data there exists 
a range of models which are favored by the data (all the 
models which lie in the degeneracy law) and so a small im- 
provement in the error bars of < Kg > (associated to the 
limited number of realizations) to better discriminate among 
the models would require a large amount of CPU time (we 
use several ALPHA DEC 3000). However, with the 4 years 
COBE data set the signal to noise ratio will improve re- 
ducing the (Q rms -ps,n) anticorrelation and therefore mak- 
ing the statistic Kg more suitable to discriminate among 
the fewer remaining models. Comparison of the Kg statistic 
with the < Xg > one shows that the former is a much better 
discriminator among different models. 



6 DISCUSSION 

The relation between Qrms—ps and n found in the previous 
section can be interpreted in terms of the coherence angle of 
the temperature random field. The interpretation is based 
on relations between mean quantities of the field which do 
not take into account all the experimental restrictions and 
the cosmic variance. Therefore this section should be only 
considered as an approximate approach to the understand- 
ing of the previous results. 

As a function of threshold level v, the mean value of 
the genus for a Gaussian random field only depends on the 
coherence angle of the field 9 C : 
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where the threshold level is given in terms of temperature 
standard deviations. The coherence angle is defined in terms 
of the ratio between the correlation function and its second 
derivative at zero lag. If it were possible to measure CMB 
anisotropies with ideal noise-less instruments, one could de- 
duce from the previous expression the intrinsic coherence an- 
gle of the underlying field (if no cosmic variance and galaxy 
cut were present). The coherence angle of noise, as inferred 
from the measured genus of the i(A-B) map and formula 
(5) is 4.3° ±0.1°. The coherence angle of the signal and noise 
53 i(A+B) map derived from its genus is 4.9° ± 0.1°. 

In the general case of a sky map including signal and 
noise 9 C is given by 
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Ci = (a?) exp(-l(l + l)a 2 eff ), 

where (a 2 ) = (& 2 m ) (see eq. 2), <J e ff is the effective Gaussian 
smoothing and Ci t s and Ci t N are the Legendre coefficients 
for the signal and noise respectively (i.e. the coefficients in 
C{9) = j^J2i( 21 + 1)CiPi(cos8)). In the extreme case of a 
noise-less map only the spectral form (i.e. n) of the model 
contributes to the coherence angle but not the amplitude of 
the signal. For instance, for the n = 1 model the DMR 
radiometers would measure an effective 9 C PS 12.4° after 
beamwidth filtering (ftj 3°), beam smearing (1.3°), smooth- 
ing angle (2.9°) and pixel size (1°, corresponding to the dis- 
persion of a Gaussian with the same area as the pixel of 
side 2.6°) have been taken into account. Beam smearing is 
the smoothing caused by the motion of the antenna dur- 
ing the 1/2 second integration time per measurement. The 
a 2 N noise coefficients can be estimated from a least squares 
fit of a harmonic expansion to the i(A-B) COBE map. A 
pure noise sky map would have a smaller coherence angle, 
9 C PS 4°. 2, considering a <J e ff given by the smoothing angle, 
pixel size and beam smearing added in quadrature. This 
value agrees with the one obtained by using equation (5). 
Moreover, equations (2) and (6) can be solved numerically 
to give the theoretical Qrms—ps as a function of n and coher- 
ence angle, Qtheovy(n; 9 C ). This law gives an anticorrelation 
between Q and n and reproduces well the empirical relation 
found in the previous section. Then, one can also estimate 
the 9 C corresponding to the 53 |-(A+B) map by fitting the 
Qtheory (n; 9 c ) to the points on the (Q rms -ps,n) plane with 
minimum < Xg >• This can be considered the coherence an- 
gle of the original random field and its value is 5.35° ± 0.1°. 
The anticorrelation found between Qrms—ps and n results 
from the fact that models with high n tend to give more 
weight to small angular scales, so in order to keep the same 
coherence angle for a given instrument noise the amplitude 
Qrms—ps of the spectrum must be decreased as n increases. 

Notice that by using equation (5) a different value of 
the coherence angle was found, as could be expected due 
to the uncertainty produced by the cosmic variance and the 
galactic cut which were not considered. The correction to 
this value can be calculated by Monte Carlo runs of a ran- 
dom field with a given 9 C = 5°. 35, and then comparing with 
the one obtained from the genus relation (eq. 5). The correc- 
tion found is just the difference between the two estimated 
values of 9 C : the value for the random field of 5°. 35 and 9 C 
estimated by eq. (5) for the COBE data. 

Finally, we point out that a similar reasoning would 
apply to the results obtained from the analysis based on the 
number of spots N , simply by substituting equation (5) by 
the following equation: 
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7 CONCLUSIONS 

Based on a generalized % 2 statistic which is defined in terms 
of the genus (or spot number) that would be measured 
by different cosmic observers we are able to constrain the 
quadrupole moment Qrms—ps and the spectral index n of 
the density fluctuation power spectrum at recombination 
with the COBE-BMR first year data. We find that the two 



Genus and spot density in the COBE DMR first year anisotropy maps 5 



parameters should lie within a region around the straight 
line Q rrns -ps = 22.2 ± 1.7 - (4.7 ± 1.3) x n for fixed n, 
which corresponds to the line of constant 9 C = 5°. 35. This 
relation has been obtained with the mean value < % 2 > and 
confirmed with the alternative quantity < K > derived from 
the Kolmogorov-Smirnov statistic applied to the \ 2 distri- 
butions. This anticorrelation is consistent with the results of 
Smoot et al. (1994), using the mean genus of the models for 
several smoothings of the data, and with the one obtained 
by Seljak and Berstchinger (1993) based on an analysis of 
the COBE correlation function. 

We have also studied the minimum uncertainty due to 
cosmic variance with which one can obtain the spectral index 
n when using the genus as the statistical quantity in the 
comparison with the COBE data. The result is a la error 
bar of Sn ~ 0.2 (for n ~ 1). 
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